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Abstract 


It is difficult to measure the WW-fusion Higgs production process (ete~ — výh) 
at a lepton collider with a center of mass energy of 240-250 GeV due to its small 
rate and the large background from the Higgsstrahlung process with an invisible Z 
(ete + hZ, Z > vv). We construct a modified recoil mass variable, m? coip 
defined using only the 3-momentum of the reconstructed Higgs particle, and show 
that it can better separate the WW-fusion and Higgsstrahlung events than the orig- 
inal recoil mass variable Mrecoil. Consequently, the MP att variable can be used to 
improve the overall precisions of the extracted Higgs couplings, in both the con- 
ventional framework and the effective-field-theory framework. We also explore the 
application of the m? coi] variable in the inclusive cross section measurements of the 


Higgsstrahlung process, while a quantitive analysis is left for future studies. 
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1 Introduction 


A lepton collider running at a center of mass energy of around 240 to 250 GeV is ideal 
for studying the properties of the Higgs boson. It can collect a large amount of Hig- 
gsstrahlung (ete~ — hZ) events, which has a cross section maximized at around 250 GeV. 
At higher energies, the WW-fusion process of Higgs production (ete~ —> vivh) can be 
better measured as its cross section increases with energy. It is important to have good 
measurements of the WW-fusion process, which provides information complementary to 
the one from the Higgsstrahlung process. In the conventional kappa framework, the WW- 
fusion process can constrain the hWW coupling and is also an important input for the 
determination of the Higgs total width.’ In the effective-field-theory (EFT) framework, 
the WW-fusion and the Higgsstrahlung processes probe different combinations of EFT 
parameters. The inclusion of both processes, as well as the diboson one (ete + WW), 
is crucial for discriminating different EFT parameters and obtaining robust constraints 
on all of them [2-4]. However, it is not guaranteed that the runs at energies higher than 
240-250 GeV will be available. The proposed Circular Electron Positron Collider (CEPC) 
in China does not have plans for the 350 GeV run at the current moment [5]. For the Fu- 
ture Circular Collider (FCC)-ee at CERN [6] and the International Linear Collider (ILC) 
in Japan [7], a significant amount of time may also be spent at the 240 GeV/250 GeV run 
before moving on to higher energies. The measurements of the WW-fusion process at 
240-250 GeV is therefore of great relevance to the study of Higgs physics. 


‘See e.g. Ref. [1] for a recent study under this framework. 


Figure 1: The two processes that contributes to the ete~ — vvh channel. Left: The 
WW-fusion process of Higgs production. Right: The Higgsstrahlung process with Z 
decaying to a pair of neutrinos. 


It is difficult to measure the WW-fusion process at 240-250 GeV for the following two 
reasons. First, it has a small rate at lower energy, with a cross section of 6.72 fb at 
250 GeV assuming unpolarized beams, while the total cross section of the Higgsstrahlung 
process is 212fb at the same energy [5]. Second, the Higgsstrahlung process with Z 
decaying invisibly (ete~ + hZ, Z — vi) is the dominate background of WW fusion. 
Both contributing to the channel ete~ — voh as shown in Fig. 1, the cross section of the 
former is more than six times of the latter one. With longitudinal beam polarizations, 
the situation is slightly better. If the signs of the polarizations are sgn(P(e~), P(et)) = 
(—,+), the WW-fusion cross section is enhanced, and by a larger factor than the one of 
ete — hZ. The method of recoil mass can also be used to separate the WW-fusion 
and hZ events, as the reconstructed mass of the neutrino pair should center around the 
Z mass for hZ events. However, the discriminating power is limited by the detector 
resolutions (especially for hadronic Higgs decays) and other effects [5, 8, 9]. A consistent 
treatment of the hZZ and hWW couplings is also required since they are related by 
gauge invariance. In the EFT framework, the relation is complicated by the inclusion 
of Dimension-6 operators which generate anomalous couplings with Lorentz structures 
different from the standard model (SM) ones. 


In this paper, we try to address the issues mentioned above and further optimize the 
measurements in ete” — vvh at 240-250 GeV. We first perform a collider study in Sec- 
tion 2 with a comparison of the recoil mass variable and its variations. We try to validate 
our study by following closely Ref. [9]. We point out that the variable moin, defined using 
only the 3-momentum of the reconstructed Higgs particle, could provide a discriminating 
power better than the original recoil mass variable Meco does. We then implement the 
Mrecoi] and M?co distributions in the EFT global analysis in Section 3, using the frame- 
work in Ref. [3]. We point out the importance of fitting the EFT parameters directly 
to the binned Myeco Or Meco distribution instead of fitting them to the extracted cross 
sections of the WW-fusion process. We also apply m?,..,, to the inclusive ete” > hZ 
process in Section 4 and comment on its potential use in the inclusive hZ cross section 
measurements. Finally, we conclude in Section 5. We provide a short summary of the 
EFT framework used in our analysis in Appendix A and the numerical expressions of the 


EFT dependence of the (modified) recoil mass distributions in Appendix B. 


The following collider scenarios are considered in our study: 


e CEPC with 5ab~' data collected at 240 GeV, with unpolarized beams [5]. This 
scenario can also be thought as the earlier stage of the FCC-ee, which also plans to 
collect 5ab~' at 240 GeV and eventually 1.5 ab! data at 350 GeV as well [6]. 


e ILC with 2ab ' data collected at 250 GeV and beam polarizations of P(e~,e*+) = 
(40.8, +0.3), which could be considered as the first stage of a full program with 
center of mass energies up to 500 GeV [10]. 


For the WW-fusion measurements, we focus on the channel with the Higgs decaying to a 
pair of bottom quarks (e+e~ + vivh,h — bb), which has the largest branching ratio. The 
measurements of WW fusion at 240-250 GeV with other Higgs decay channels are not 
reported in the official documents due to the poor constraints (see e.g., Refs. [5, 10, 11)). 


2 The modified recoil mass of ete~ > vih,h — bb 


At lepton colliders, the method of recoil mass can be used to reconstruct the mass of a 
particle without measuring its decay products. One of its most important applications 
is the measurement of the inclusive rate of the Higgsstrahlung process, ete — hZ. 
Assuming both the Higgs and Z are on mass shell, one could write the relation 


mp = Ep — |pn|? = (Vs — Ez)? — pz’, (2.1) 


where the total center of mass energy ys is fixed up to corrections from beam energy 
spread and initial state radiations. By measuring the energies and momenta of the Z 
decaying products one could reconstruct the mass of the Higgs particle. This can be 
used to select ete” —> hZ signal events without tagging the Higgs decay products, which 
makes it possible to measure the inclusive cross section of this channel. It also provides 
the best Higgs mass measurement. For example, a precision of 5.9 MeV can be achieved 
with the leptonic Z decay channels of the inclusive hZ measurements at the CEPC [5]. 


If the Higgs decay products are measured, the recoil mass can be turned around to 
reconstruct the Z mass, since the following relation also holds for an ete” — hZ event, 


mz = (Vs — En)” — |pal? = s — 20/8 En + mp. (2.2) 
The recoil mass can then be defined as 


Me scoi] = S — 2/8 Er + (mi)? , (2.3) 


recoil 


rec 


where E} and m}® are the reconstructed Higgs energy and mass. For Higgs decaying 
to a pair of bottom quarks, they are the total energy and the invariant mass of the 
two b-jets. This offers a way to separate the Higgsstrahlung events with an invisible Z 
(ete + hZ, Z > vv) from the WW-fusion events, both contributing to the channel 
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ete” — vvh, as shown in Fig.1. However, due to finite jet resolutions, beam energy 


spread and other effects, the recoil mass distribution of the hZ events has a rather large 
spread. This limits it discriminating power especially at the energy 240-250 GeV, for 
which the recoil mass distribution of the WW-fusion events spreads around the same 
region. We also find that for Higgs decaying to a pair of b-jets (h — bb) the uncertainty 
on the recoil mass is dominated by the energy and momentum resolutions of the b-jets. 
This is also obvious from the observation that the recoil mass distribution for the Higgs 
mass reconstruction in Eq. (2.1) is much narrower for the leptonic Z decay channel than 
for the hadronic one (see e.g. Ref. [5]). 


While the recoil mass defined in Eq. (2.3) makes use of all the kinematic information, 
it does not make any assumption on the value of the Higgs mass. Both the Higgs width 
and the projected uncertainty of its mass are at the MeV level and can be neglected 
compared with the effects of jet resolution. Using the information of the Higgs mass, two 
modifications of the recoil mass can be constructed. The first, using only the reconstructed 
Higgs energy as the measurement input, is defined as 


Gites) = §— 2/8 Ex + mi, , (2.4) 


where mp is fixed to be the Higgs mass, 125.09 GeV. The other, using only the recon- 
structed 3-momentum (p;,°°) of the Higgs, is defined as 


Cee eg ae oe 2/s mi, + | Prel? + mi, ’ (2.5) 


where my, is again fixed to be 125.09 GeV. At the truth level, Mrecoi, MEcon and MPi are 
all equivalent. However, the uncertainties in the energy and momentum measurements 
certainty have different impacts on the three variables. To illustrate this impact, we 
define, for a given event, a set of five parameters {6m, 6%, ô? , Oz, Ôp} which parameterize 
the differences between the reconstructed quantities and the true ones, with 


™Mrecoil = mMysooil (1 + ô m), it ce = Merseoil (1 + Sn) > Meecoil = Miil T OF a) ’ (2.6) 
and 

Ex’ = E,(1+ de), Pa | = [Pall + op), (2.7) 

where mitt. is the true parton level recoil mass, Ep and pp, are the true energy and 3- 


momentum of the Higgs. For hZ events, mt% = mz (assuming it is on shell), and the 


three parameters ôm, 6% and 6?, can be written in terms of dg and 6,. At leading order, 
they are given (for hZ events) by 


1 5 

Ôm © — a (Vs — En) En ôg + RA , 
Z 
s 

ôE x — s En On, 

rae _ vs Pri? 5 


m3 En 


Sp. (2.8) 
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Figure 2: Left: The distributions of 6g and 6, (defined in Eq. (2.7)) for ete™ + hZ, Z > 
vv,h — bb at CEPC 240 GeV, after applying a Higgs-mass-window cut of 105GeV < 
mie < 135GeV. Right: The corresponding distributions of Myeco, M~.4 and Moa, 


defined in Eqs. (2.3-2.5). A radius of R = 1.5 is used in the jet clustering algorithm. The 
details on the simulations and selection cuts are stated later in this section. 


Note that dg and 46, can be either positive or negative. The overall negative coefficients 
in Eq. (2.8) indicates that if the measured energy or 3-momentum of the Higgs is larger 
than its actual value, the recoil mass variables will be smaller than the Z mass, and vice 
versa. For a fixed center of mass energy (ys = 240 GeV or 250 GeV), the values of Ep, 
and |p;,| are fixed. In particular, near the hZ threshold |p| is significantly smaller than 
En. With |p,| ~ 51GeV and Ep œ 135 GeV at ys = 240 GeV, and |p;,| ~ 62 GeV and 
E, © 140 GeV at ys = 250 GeV, Eq. (2.8) thus becomes 


in | E | 8 eee — 0.326, / — 3.9 őr / — 0.575, at 240 GeV (2.9) 
ee —1.9d¢ — 0.466, / —4.26¢ / — 0.836, at 250 GeV ’ 

where the small coefficients of 5, come from a suppression factor of ~ |pj,|?/E? relative to 
the ones of ôg, shown in Eq. (2.8). The distributions of ôp and 6, for ete™ —> hZ, Z > 
vv,h — bb at CEPC 240 GeV are shown on the left panel of Fig.2, after applying the 
selection cuts which include a Higgs-mass-window cut of 105 GeV < m} < 135 GeV on 
the b-jet pair. While 6, has a slightly larger spread than ôg, its coefficients in Eq. (2.9) are 
much smaller. We therefore expect the distribution m?...,, to have the smallest spread, 
and the one of mË oi to have the largest. This is verified on the right panel of Fig. 2 where 
the distributions of Mrecoil, MEcoi and MPi are shown. For the WW-fusion events, we 
expect a less significant difference among the distributions of the three variables (which 
are shown later in Fig. 3 & 4), since they do not have a Z in the event. The corresponding 


distributions for ILC 250 GeV are very similar to the ones in Fig. 2. 


It should be noted that the distribution of ôg in Fig. 2 is asymmetric, suggesting that 
on average the measured energy of the b-jet pair is smaller than its actual value. This is 
due to the fact that in our simulation we do not apply any jet energy corrections that are 
widely used in the LHC experiments [12, 13]. As a result, the central values of the Myecoi 
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and m£.,, distributions are also shifted to be larger than mz. Assuming a jet energy 
correction mechanism will be implemented in the future lepton collider(s), one would 
expect a corrected central value and also some improvements on the energy measurement, 
and the Meco (and m~,,,,) distribution will have a peak value around mz and a slightly 
smaller spread. We do not expect the lack of jet energy correction to have a significant 
impact on our results since the m? eoi] distribution still has a much smaller spread due to 
the parametric suppression of it uncertainty near the hZ threshold as discussed above. ” 

Having found that the variable m? coi, could better reconstruct the Z mass than Myecoil; 
we perform an analysis based on a fast simulation to explicitly exam their discriminating 
powers on the WW-fusion and hZ events at the CEPC 240 GeV (with unpolarized beams) 
and ILC 250 GeV (assuming P(e~,e+) = (—0.8,+0.3)). We generate events for both 
processes using Madgraph5 [14], which are showered with Pythia [15] before passing to 
Delphes [16] with ILD cards (using the detector geometry and flavor tagging efficiencies 
given in Ref. {17]) for detector simulations. The interference term between the WW-fusion 
and hZ processes are ignored. It should be noted that the effects of ISR photons are not 
considered in the simulation with Madgraph5. However, we expect their effects to be much 
smaller compared with the ones of jet resolutions. We use the ILC analysis in Ref. [9] 
as a guide to validate our results from the simple simulation. While the Durham jet 
clustering algorithm is used in Ref. [9], it pointed out that the anti-k, jet algorithm with 
jet radius R = 1.5 has a similar performance in the Higgs invariant mass reconstruction, 
which is used in our simulation. We also follow closely the selection cuts in Ref. [9]. In 
particular, each event is required to have exactly two b-jets and a cut on the invariant 
mass of the b-jet pair, 105 GeV < mj < 135 GeV, is applied to reduce the backgrounds. 
The cuts related to variables in the Durham jet clustering algorithm are replaced by the 
simple requirement on jet number (= 2) in each event. After event selections, we scale the 
number of signal events of ILC 250 GeV to the ones in Ref. [9] (normalized to 2ab7'). A 
similar scaling is also applied for CEPC, taking count of the differences in cross sections 
and selection efficiencies between CEPC and ILC. 


The composition of the background in the ete~ — vvh channel is also listed in 
Ref. [9]. The major ones are vibb and qq, which contributes to 42% and 34% of the 
total background after selection cuts. The gq background is difficult to simulate due to 
its huge cross section and tiny selection efficiency. For simplicity, we simulate only the 
vibb background, apply the selection cuts and scale it up to match the total background 
number, given in Ref. [9] and normalized to our run scenarios. We expect this simple 
treatment to provide a reasonable estimation on the effects of the backgrounds. 


After selection cuts, the Meco and MË oi, distributions of hZ (Z > vv), WW-fusion 
and background events are shown in Fig. 3 for ILC 250 GeV (2ab~' data with P(e~,e+) = 
(—0.8, +0.3)) and Fig. 4 for CEPC 240GeV (5ab™' data, unpolarized beams). In Fig. 3 
and Fig.4, the Mrecoi (™».,9;) distribution is shown on the left (right) panel, while the 
distributions in the bottom panels are simply the amplified versions of the ones in the 


2We thank Zhen Liu for very valuable discussions on the topic of jet energy corrections. 


ILC 250GeV (2/ab) 


1500 - — hz, 257 | 
f === WW fusion | 2500 === WW fusion | 
sanse background 4 ====- background] 


500 - 


Pig 
Pe sn a, 
= rion Aloe 


Mrecoil [GeV] 


40 60 80 100 120 140 160 


ILC 250GeV (2/ab) 


— hZ,Z>w 


(iF ent [GeV] 


ILC 250GeV (2/ab) ILC 250GeV (2/ab) 
800 - l ' 800} i 
L =— hZ, Z>w — hZ, Z>VV 
| === WW fusion === WW fusion 
600; aasia background oe 600F ===- background 
H~ 400+ papt! 
rT s ri ta 
i | i. 


200+ 


40 60 80 100 120 140 
mbecoil [GeV] 


40 60 80 100 

Mrecoil [GeV] 

Figure 3: The mrecoi (left panel) and m?..,;; (right panel) distributions of hZ (Z > vv), 

WW-fusion and background events after selection cuts for ILC 250 GeV with a luminosity 

of 2ab~' and beam polarization P(e~,e+) = (—0.8,+0.3). The distributions in the 

bottom panels are the amplified versions of the ones in the top panels (and also with a 
different bin size). 


top panel. To check the validity of our results, we compared our Meco) distributions to 
the ones in Refs. [9, 5] and found a reasonable agreement in terms of spread ranges and 
distribution shapes. Comparing the mreco and M?co distributions in Fig.3 and 4, it is 
clear that m?..;, provides a better discrimination between hZ and WW-fusion events, 
with the m? coi distribution of hZ having a much sharper peak around the Z mass than 
the Mrecoi) One. On the other hand, the m” coi distribution of the vvbb background has a 
more even spread than the Mrecoi one. This is because that for the vibb background the 
bb pair does not come from the Higgs decay. Using the wrong mass assumption therefore 
makes the reconstruction of Z mass worse. 


Following Ref. [9], we apply a x? fit to the binned myecoi, and MP ooi distributions to 


extract the precisions (one-sigma uncertainties) of the hZ and WW cross sections, denoted 
as Ong and owwa. This is done by treating the overall rates, opz and oww n as free 
parameters while assuming a perfect knowledge on the the shapes of the distributions. 
Ref. [9] also treats the overall cross section of the background (opg) as a free parameter. 
We consider two cases with op, treated as a free parameter as well as fixing it to the 
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Figure 4: Same as Fig. 3 but for CEPC 240 GeV with 5ab~! data and unpolarized beams. 


P 


predicted value. The total y? of the Mrecoil or MËoi distribution is given by 


x? — 5 (ico ~~ He) 


i 
a Nexp 


l (2.10) 


where for each bin 7, i is the expected number of events from simulation and Pos is a 
function of onz and oww (and opg). To ensure enough statistics in each bin, we choose a 
bin size of 5 GeV except for the first and last bin, which are chosen to include all the events 
below 65 GeV (60 GeV) and above 130 GeV (115 GeV) for the Mrecoil (™?,¢;) distribution 
at ILC, and all the events below 65 GeV (60 GeV) and above 120 GeV (110 GeV) for the 
Mrecoil (Maoi) distribution at CEPC. The results of the x? fits are presented in Table 1 for 


ILC 250 GeV (2ab7' data with P(e~,e+) = (—0.8, +0.3)) and Table 2 for CEPC 240 GeV 
(Sab ' data, unpolarized beams). 


Overall, our results on the precision of the WW-fusion cross section is slightly worse 
than the ones in Ref. [9] (if normalized to the same luminosity) and the CEPC preCDR [5]. 
This is not surprising, since the results could depend on details of the simulation, for which 
we only performed a simplified study. In what follows, we shall focus on the relative 
difference between the results from Mreco and Meco While the distributions of meoil 


clearly better separates the hZ and WW-fusion events as shown in Fig.3 and Fig. 4, 
in the 3-parameter fit (shown on the left panels of Table 1 and Table 2) the precision 
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ILC 250 GeV, 2ab~‘, uncertainties normalized to SM predictions 
3-parameter fit fixing Op, 
Mrecoil uncertainty correlation matrix uncertainty | correlation matrix 
OhZ OWW=>h Obg OhZ OWW—>h 
ORZ 0.049 1 0.47 -0.97 0.011 1 -0.69 
OwW-—h 0.063 1 -0.63 0.045 1 
Tbg 0.31 1 
3-parameter fit fixing Obg 
MP coil uncertainty correlation matrix uncertainty | correlation matrix 
OhZ OWW=>h Obg ORZ OWW—>h 
ORZ 0.010 1 0.21 -0.51 0.0088 1 -0.46 
OWW—>h 0.059 1 -0.83 0.033 1 
Obg 0.088 1 


Table 1: The one sigma uncertainties and correlations of the cross sections of the Hig- 
gsstrahlung process with an invisible Z (onz) and the WW-fusion process (oww_sp) at 
ILC 250 GeV from a fit to the Mrecoi (top panel) and Moi (bottom panel) distribu- 
tions. A total luminosity of 2ab~' with beam polarization of P(e~,e+) = (—0.8, +0.3) 
is assumed. In the 3-parameter fit on the left panel, the overall normalization of the 
background is treated as a free parameter in the fit. In the 2-parameter fit on the left 


panel, the total number of background events is fixed to the predicted value. 


of owwa from the m?.., distribution is similar to (or even worse than) the one of 
Mrecoil, due to the fact that WW-fusion and background events have a larger overlap in 
the moi distribution than in the Mrecoi one. Meo nevertheless significantly improves 
the precisions of opz and opg in the 3-parameter fit. Assuming a good knowledge of 
the background, one may also fix the background cross section to the predicted value. 
In the 2-parameter fit with opz and oww-sn, the Meo distribution indeed provide a 
significantly better constraint on owwa, with an improvement of about 30% at ILC and 
20% at CEPC, compared with the constraint from the Myecoi distribution. 

In the kappa framework with the hZZ and hWW couplings treated as independent 
parameters, the WW-fusion measurement is an important input for constraining the 
hWW coupling (with the other being the Higgs decay, h > WW*). With a 20-30% 
improvement on the precision of the WW-fusion cross section, a sizable improvement for 
the constraint on the hW W coupling is expected. The WW-fusion cross section is also an 
important input for the determination of the Higgs total width, following the relation [5]? 


x T(h > bb) m a(vih, h — bb) 
BR(h > bb)  BR(h > bb) - BR(h > WW*) ’ 


3 CAEN 


where, following the usual convention, o(vvh) denotes only the WW-fusion contribution 


3It should be noted that, despite the usual claim of being model independent, Eq. (2.11) explicitly 
assumes that the hWW coupling is independent of the energy scale (i.e., anomalous couplings such as 
hW"’W,, are absent), which is not true under the more general EFT framework. 
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CEPC 240 GeV, 5ab~‘, uncertainties normalized to SM predictions 
3-parameter fit fixing Op, 
Mrecoil uncertainty correlation matrix uncertainty | correlation matrix 
OhZ OWW=>h Obg OhZ OWW—>h 
ORZ 0.024 1 0.28 -0.95 0.0077 1 -0.61 
OWW-— sh 0.058 1 -0.47 0.051 1 
Tbg 0.20 1 
3-parameter fit fixing Obg 
MP coil uncertainty correlation matrix uncertainty | correlation matrix 
OhZ OWW=>h Obg ORZ OWW—>h 
OhZ 0.0071 1 0.098 -0.35 0.0066 1 -0.45 
OwWw-— sh 0.083 1 -0.87 0.041 1 
Obg 0.082 1 


Table 2: Same as Table 1 but for CEPC 240 GeV with 5 ab‘ data and unpolarized beams. 


to it. With a 20-30% improvement on the precision of o(vvh, h — bb), as well as some 
possible improvement on the determination of BR(h — bb) from a better measurement. of 
a(hZ, Z = vv, h = bb), the precision of the Higgs total width obtained using Eq. (2.11) 
could be improved by at least 20-30% using the m? coi] variable. 


To conclude this section, we would like to emphasize that, while we try our best to 
validate our results, they do rely on simple simulations and should be explicitly tested by 
experimental groups with proper simulation tools. The fits performed in obtaining the 
results in Table 1 and 2 also assume a perfect knowledge of the distribution shapes for 
each process, which may not be a good assumption in an actual experiment. We also 
include only the vibb background, while other backgrounds may have different kinematic 
features. Nevertheless, we expect our results to still hold qualitatively due to the simple 
reasoning that m? coi has a smaller uncertainty and better reconstruct the Z mass. 


3 Improving Higgs coupling constraints in the EFT 


framework 

Having explored the capability of the m? oi variable in improving the measurement of 
the WW-fusion process, we are now ready to exam its impact on the determination of the 
Higgs couplings. We choose to study the Higgs coupling constraints in a global effective- 
field-theory (EFT) framework with dimension six (D6) operators. Such a framework 
has several advantages. First, assuming the scale of new physics is high, the EFT with 
D6 operators gives a good parameterization of the effects of new physics and the results 
can be mapped to any specific model that satisfies the assumptions of the framework. 
Second, it takes count of the connections among different measurements. For instance, 


“For recent Higgs EFT studies in the contexts of future lepton colliders, see Refs. [2-4, 18-23]. 
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some operators contribute to both Higgs processes and the dibson process, and the triple 
gauge coupling (TGC) measurements from the diboson process can thus help the overall 
constraints on the Higgs couplings [24]. Gauge invariance is also imposed by construction 
in the EFT framework. We focus on the CEPC 240GeV and follow Ref. [3] in terms 
of the basis choice and measurement inputs. In particular, focusing on the Higgs and 
diboson measurements at 240 GeV, and making reasonable assumptions, a total number 
of 11 parameters are sufficient to describe the contributions from the D6 operators. A 
short summary of the framework in Ref. [3] is provided in Appendix A. The methods we 
propose should nevertheless be applicable to other collider scenarios and frameworks. 


A few important differences between the cross-section fit in Section 2 and the EFT 
analysis should be noted. While not specifically mentioned, the cross-section fit does make 
assumptions on the new physics, in particular that it only modifies the overall rates, not 
the differential distributions, of the hZ and WW-fusion processes. In the cross-section 
fit, the hWW and hZZ couplings are also assumed to be independent, regardless of 
their relation from gauge invariance. The EFT analysis, while imposing gauge invariance, 
contain anomalous couplings of the form hZ,,,Z'” and hZ,0,Z"” (and the same for W) 
which have different momentum dependences from the SM couplings due to the extra 
derivatives. The potential new physics contribution to the hZy vertex is also included in 
the EFT analysis, which could contribute to ete” —> hZ via an s-channel photon. It is 
thus beter to directly fit the EFT parameters to the Myecoii or M?,.,, distribution of the 
inclusive ete” — vvh process instead of fitting them to the extracted precisions of cross 
sections from Section 2. By fitting to the inclusive ete~ — výh process we also include 
the interference term of WW-fusion and hZ processes, which is usually ignored in cross- 
section fits. The expressions for the total cross section of ete~ + vvh and the (binned) 
differential ones of the Myecoi and MP eoi] distributions in terms of the EFT parameters 
are listed in Appendix B. 


The measurement inputs of the Higgsstrahlung (ete~ — hZ) and diboson (ete~ > 
WW) processes at CEPC 240 GeV are listed in Table 3. The estimations of hZ mea- 
surements are taken from Ref. [25], which updates the ones in the CEPC preCDR [5]. 
In addition, The angular observables of ete~ — hZ in Ref. [19] are included, for which 
we use only the channel ete~ > hZ, h > bb, Z — C+ and assume a fixed 60% signal 
selection efficiency, following Refs. [21, 3]. For the TGC measurements, we follow the 
treatment in Ref. [3] which adopts the one in Ref. [26] with the addition of a universal 1% 
systematical uncertainty in each bin of all differential distributions. We directly list the 
resultant one-sigma constraints of the anomalous TGC parameters and their correlations 
in Table 3. We construct the total y? by summing over the 7s of all measurements and 
perform global fits to obtain the precision reaches (one-sigma bounds) of the relevant EFT 
parameters. 


We consider three scenarios in the global analysis. All three use the inputs on Hig- 
gsstrahlung and TGC measurements in Table 3, but different information on the measure- 
ment of ete” — výh. The first one uses only the total rate of ete — výh. The second 


(third) uses the information in the Myecoi (Maoil) distribution, with the EFT parameters 
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CEPC 240 GeV, 5ab~‘, unpolarized beams 
ete > hZ ete > WW 
olete —> hZ) 0.50% uncertainty | correlation matrix 
a(hZ) x BR 691.2 OKy AZ 
h — bb 0.24%* 6g1,z | 6.4 x 107° 1 0.068 -0.93 
h> c 2.5% k, | 3.5 x 1073 1 -0.40 
h > gg 1.2% Az | 6.3 x107? 1 
h> tr 1.0% 
h => WW* 1.0% 
h —> ZZ* 4.3% Angular observables in 
h> yy 9.0% ete” > hZ, h > bb, Z > He 
h> pp 12% are also included. 
h> Zy 25% 


Table 3: A summary of the measurement inputs from ete~ + hZ and ete” + WW at 
CEPC 240 GeV used in the EFT fit, assuming a total luminosity of 5 ab~' and unpolarized 
beams. Inputs on rate measurements of e*e~ — hZ are from Ref. [25] which updates the 
estimations in the preCDR [5]. For the precision of o(hZ) x BR(h — bb) (marked by a 
star *), we have excluded the contribution from ete~ > hZ,Z > vv,h > bb to avoid 
double counting with ete~ > vih,h > bb. The angular observables in ete~ > hZ, h > 
bb, Z — 4+6 are included, assuming a 60% signal selection efficiency. The constraints 
on aTGC parameters from measurements of ete” — WW are obtained following the 
treatments in Ref. [3]. 


directly fitted to the binned distributions. We also compare the reach with the more 
conventional method of fitting the EFT parameters to the extracted precisions of hZ and 
WW-fusion cross sections in Section 2, ignoring the correlation between the two cross 
sections (which is often not reported). Since the background is assumed to be SM-like in 
the EFT analysis,’ for the extracted precision of cross sections we also use the results of 
the two-parameter fit with fixed background, shown on the right panel of Table 2. The 
results of the 11-parameter fit are presented in Fig. 5. 


Comparing the reaches of the three scenarios (shown by the orange, blue and cyan 
columns in Fig.5), we indeed observe a non-negligible improvement on the overall reach 
if the information in the Mrecoi Or ™Me..o; distribution is used. In particular, the reach 
on the parameter dcz (corresponding to a shift in the SM hZZ and hWW couplings) is 
improved by more than 10%. Consequently, the reach on ce OYc, OY» and dy,, which 
contribute to the Higgs decay to gg, cé, bb and 77, have also been improved by a similar 


(or slightly less) factor. The reach with the m?,.,,, distribution is better than the one with 


5Tt is reasonable to fix the background (which has no Higgs) to the SM predictions in an EFT global 
framework in this case, as deviations from SM are strongly constrained by the electroweak precision 
measurements at Z-pole or other measurements. Fixing the background nevertheless requires one to 
have a very good knowledge of the total rate and distribution shape of the background as pointed out in 
Section 2. 
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Figure 5: The one-sigma precision reach of the 11-parameter fit in the EFT framework 
at CEPC 240GeV with 5ab~' data and unpolarized beams. See Appendix A for the 
definitions of parameters. Three scenarios are shown which differ on the information used 
for the ete” — výh measurement. The first one uses only the total rate of ete > voh. 
The second (third) one uses the Myecoit (Mco) distribution, and contain two sub-scenarios. 
The one shown by the light shaded columns uses the extracted precisions of opz and 
Oww-— nr in Section 2 (the two-parameter fit on the right panel of Table 2, correlation 
ignored). The one shown by the solid columns is obtained from a direct fit to the binned 


Mrecoil (M9) distribution. 


Mrecoil aS We expected. However, the relative improvement from Myecoi to Moco turns out 
to be very marginal. We also find that, if fitting the EFT parameters to the extracted 
cross sections caz and oww-s, without taking count of their correlation (the results are 
shown with light shades for the 2nd and 3rd columns), the reaches are worse than the 
ones from direct fitting the EFT parameters to the distributions, in particular for the 
Mrecoil distribution. This is because, as shown in Section 2, the uncertainties of opz and 
Oww- >n have a large correlation between them due to the difficulty in separating the two, 
in particular for the Myecoy distrubiton. This correlation is usually not reported in official 
documents, and the omission of it could lead to a considerable impact on the overall 
reach. It should be noted that, our results using only the total rate of ete — vih is 
worse than the corresponding ones in Ref. [3]. This is because our estimations on the rate 
measurement of ete~ — vvh is more conservative than the one in Ref. [3], which is derived 
from the CEPC preCDR [5]. If the overall cross section measurement of ete% — voh 
can be improved (e.g. by optimizing the selection cuts), we expect the use of Mrecoil 
and m?..,,,; distributions would also bring a more significant improvement on the overall 
reach of the EFT fit. We have also chosen very conservative bin sizes to control the 
uncertainties in each bin from simulation. Further optimizations of the analysis may also 


provide substantial improvements on the reach with the myecoy and MË coi distributions. 
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We also find that, if the TGCs can be measured with much better precisions such as 
the ones in Ref. [4] for ILC 250 GeV (which are one order of magnitude better than the 
ones in Table 3), or if multiple runs with different beam polarizations are available (also 
likely to be the case for ILC), the improvement from the Mmyecoy and Moco distributions 
with respect to using only the total rate of ete — vvh becomes rather insignificant. 
This is not surprising, since very precise TGC measurements can effectively remove two 
degrees of freedom in the fit so that there is less need for additional handles to discriminate 
the parameters. The interference term of the ete” + hZ diagram with an s-channel Z 
and the one with an s-channel photon is also sensitive to the beam polarization, which 
can help probe the operators that contribute to this interference [3, 4]. In general, once 
a sufficient number of constraints are included in a global analysis, the overall precision 
reach is expected to be less sensitive to the impact of a single measurement, such as the 
one of ete” — vvh. It is nevertheless important to optimize the measurements in order 
to maximize the sensitivity to new physics. 


4 Applications on the inclusive hZ measurements 


While we have focused on the WW-fusion measurements in the previous two sections, a 
question of great interest is whether the variable m? soi could be applied to improve the 
inclusive measurement of ete~ — hZ, where the decay product of Z are tagged instead. 
While suffering from the jet resolution, the hadronic Z decay channel provides a slightly 
better measurement on o(Zh) than the leptonic one, thanks to its large branching ratio.° 
An improved measurement of the hadronic Z channel could thus have a significant impact 
on the overall precision reach of o(hZ). It is straightforward to write down the recoil mass 
and its two variations for the reconstruction of the Higgs mass in ete~ —> hZ, which are 


ieil TS 2V5 Ex + +(m aT, 


Cans = 3 = 2\/s Ez a me; 
Caen = g= 2y/s mz + | pr? + my ’ (4.1) 


where E'7°, pz° and m° are the reconstructed energy, 3-momentum and invariant mass 
of the Z ed the two jets, while mz is the true Z mass, fixed to be 91.19 GeV. Similar 
to Eq. (2.9), we derive the deviations in the measured Myecoi, Me oo; Maoi aS a function 
of the deviations in the measured Z energy and 3-momentum to be 


fn J oE / a m {O91 bn 0.176, / -168g / —0.396, ab MOG, 
mi ml °m™~) 0996 —0.255,/ —186n / —0.566, at 250Gev? 
where ôm, 6& and ô, are defined as 
Mrecoil = Mrecoilll + ôm),  Miecoil = Mees + ôn) Mecon = Mrecsa(l + ôR), (4.3) 


°For instance, the precision of the inclusive hZ cross section measured from the leptonic (hadronic) 
Z channel is reported to be 0.8% (0.65%) in the CEPC preCDR [5]. 
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Figure 6: The distributions of ôg and 6, of the reconstructed Z (defined in Eq. (4.4)) in 
ete — hZ, Z — qq at CEPC 240 GeV after applying a Z-mass-window cut of 70 GeV < 
Mqq < 110 GeV (left) or 80GeV < Mag < 95 GeV (right). The Higgs is forced to decay 
invisibly in the simulation to ensure the correct reconstruction of Z. A radius of R = 0.5 
is used in the jet clustering algorithm. 


with mit’. = mp. For dg and 6,, the definitions are 
E% = Ez(1+ dz), Pz" | = |Pal(1 + 4p), (4.4) 


where Fz and pz are the true energy and 3-momentum of the Z. Similar to Eq. (2.9), in 
Eq. (4.2) the coefficients of 6, is also smaller than the ones of dg, but with a suppression 
factor of ~ |pz|?/Ez instead. The distributions of ôg and 6, for the reconstructed Z are 
shown in Fig. 6, with the details of simulation stated later in this section. Note that the 
cut on the Z-mass window has a strong impact on the distributions of dg and ôp. For 
larger deviations of the measured energy and momentum from the true ones, the invariant 
mass also tends to be further away from its true value. 


To compare the reconstruction power of Mrecoi and m?,..., on the Higgs mass, we 
perform a simple analysis using the simulation tools listed in Section 2. One important 
difference here is that for Higgs inclusive measurement with Z — qq, the final states 
could contain additional jets from Higgs decay, making it more difficult to reconstruct 
the Z. Due to the additional jets, we set the jet radius to R = 0.5 in order to reduce 
the contamination among the jets. For an event with more than two jets, we choose 
the pair of jets with an invariant mass that is closest to the value of Z mass. We then 
apply a Z-mass-window cut on the invariant mass of the jet pair, mg, intended for 
removing backgrounds. The difference between mMreco and m»..;, is strongly correlated 
with the size of the Z window — in the limit that the invariant mass equals the actual Z 
Mass, Mrecoi) and m?...;, become equivalent. We therefore consider both a larger window, 
70 GeV < Mag < 110 GeV, and a smaller one, 80 GeV < mg < 95 GeV. The distributions 
of Mrecoit and M”,..;, for ere” + hZ, Z — qq after the selection cuts are shown in Fig. 7 for 
CEPC 240 GeV. To estimate the impact of the combinatorial problem in the reconstruction 
of Z, we first consider a case in which the Higgs is forced to decay invisibly in the 
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Figure 7: Top row: The distributions of Mrecoi and m?..,,, for ete” > hZ,Z > qq 
at CEPC 240GeV. The Higgs is forced to decay invisibly in the simulation to avoid 
the combinatorial problem. Bottom row: The same distributions with Higgs inclusive 
decay. The left and right panels differ on the Z-mass-window cut, which is 70GeV < 


Mgqq < 110 GeV (80 GeV < Mag < 95 GeV) for the left (right) panels. 


simulation. The only purpose of the invisible decay is to avoid having additional jets 
from the Higgs decay and ensure a clear identification of the Z jet-pair. The results are 
shown on the top panels of Fig. 7 for the two choices of Z-mass-window cuts. For this ideal 
case, it is clearly that m?o has a significantly narrower spread and provides a much better 
reconstruction of the Higgs mass than Mrecoi does. The improvement with Mm? coi] is more 
significant if a large Z-mass window cut is applied as we expected. For the realistic case 
with Higgs inclusive decays, the distributions are shown in the bottom panels of Fig. 7. 
The reconstruction of the Higgs mass is worse for both myecoi and m?,.,,, distributions 
due to the wrong jet-pairing. However, m?.,,;, still has a better performance than myecoil, 
so its usefulness is not washed out by the combinatorial problem. We also note that, due 
to the lack of jet energy correction mentioned in Section 2, our distributions of Myecoil 
peak around 130 GeV rather than 125 GeV. While the central values of the distributions 
can be corrected, we expect Maoi to still have a better performance than Mrecoi after 
the implementation of jet energy corrections due to the parametric suppression on the 
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Figure 8: The distributions of reco and m?,.., for WW- —> qqqg at CEPC 240 GeV 
with a Z-mass-window cut of 70GeV < mg < 110 GeV (left panel) or 80 GeV < Mag < 
95 GeV (right panel). 


uncertainties of M?co near the hZ threshold. 


Since the background events do not have Higgs in them, we do not expect the m?,..., 
distribution of them to accumulate around the Higgs mass. As a simple estimation, we 
show the reco and MË coi) distributions for one of the main backgrounds, WW- —> qqqq 
in Fig 8, also for both choices of the Z-mass-window cuts. It is interesting to notice that 
for the larger Z mass window, m?._,;, actually reduces the paces round events in the region 
of ~100-120 GeV, while for the smaller window, the Mrecoi and mË distributions are 
very similar. 


recoil 


Our study shows that the m?o variable could better reconstruct the Higgs mass 
for the signal and does not have the same effect on backgrounds. As such, we expect 
it to provide a significant improvement on the inclusive cross section measurements of 
the Higgsstrahlung process compared with the conventional recoil mass variable Myecoil- 
Needless to say, such an improvement is crucially relevant to the studies of the Higgs 
boson properties. We also find a similar behavior of the signal and background distri- 
butions at the ILC 250 GeV, the results of which are not specifically shown. Since we 
have only performed a simplified simulation analysis and have not considered some of the 
important backgrounds, we will restrain ourselves from doing any quantitative analysis 
on the inclusive 0(hZ) measurements and leave it for the experimental groups who have 
better tools for such an analysis. 


5 Conclusions 


In this paper, we explore the use of the recoil mass and its variations in the measurements 
of the WW-fusion procese at a lepton collider with a center of mass energy of 240-250 GeV. 
We found the variable m?.,...,, constructed using only the 3-momenta of the Higgs decay 
products, can better separate Higgsstrahlung events with an invisible Z from the WW- 
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fusion events than the original recoil mass Meco does, with an improvement up to 20- 
30% on the precision of the WW-fusion cross section. We study its impact in both the 
conventional framework and the effective-field-theory one. In the conventional framework, 
a better precision on the WW-fusion cross section leads to a significant improvement on 
the constraints of the hWW coupling and the total Higgs width. In a global analysis 
under the effective-field-theory framework, using the information in the Meco Or M?,.531 
distributions could improve the reach on some of the EFT parameters by more than 
10% compared with just using the total rate of the ete — vvh channel. We find that 
fitting the EFT parameters directly to the binned distributions gives the best precision 
reach. On the other hand, if the EFT parameters are fitted to the precisions of the WW- 
fusion and AZ, Z — vv cross sections extracted from the Mmyecoy distribution, the precision 
reach could suffer from the large correlation between the two cross sections if it is not 
taken count of. We also explore the use of m?.,.., in the inclusive measurements of the 
Higgsstrahlung process (ete~ — hZ) with hadronic Zs and find that it can significantly 
improve the reconstruction of the Higgs at a center of mass energy of 240-250 GeV. The 
use of M? coi] could therefore potentially lead to an improvement on the overall precision of 
the inclusive hZ cross section measurements. The construction of moi is also extremely 
simple and does not require any additional measurements. It should be straightforward 
to implement mM? in any studies that make use of the recoil mass distribution. 
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A The “12 (or 11)-parameter” effective-field-theory 
framework 


We follow the framework in Ref. [3] which uses the Higgs basis, proposed in Ref. [27] 
and applied also in the studies of LHC Higgs measurements in Refs. [28, 24]. We focus 
on CP-even dimension-6 (D6) operators and omit the ones that induce fermion dipole 
interactions. We also assume the Z-pole observables and W mass to be SM-like, given 
that they are already very well constrained by LEP and can be further constrained with 
a Z-pole run at the future lepton colliders. 


The relevant parts in the Lagrangian of the SM and D6 operators are 


LD Lavy + Laff + Lige , (A.1) 
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where the Higgs boson couplings to a pair of SM gauge bosons are given by 


2-2 2 12\2:2, 
a WIW + (1+ 6cz) (g 1 jë 


h 
Lavy = a (1 + dcw) Zu Z" 


2 
+ eww E WEW + cwo g? (Wr AW! + hc.) 
2 2 [72 fh 
s na a uv e v e + 7 
T Cgg a k F Coy qa At F CZy g g Z pA" 
2 12 


T ZZ E ya + cz g Z0 Z” + Cy 99 Z ð, A” : (A.2) 


The parameters in Eq. (A.2) are not all independent. Four constraints can be written 
down by imposing gauge invariances, which we choose to rewrite dcw, Cww, Cwo and cy 
as 


dcw = 6cz + 46m, 


cww = Czz + 255, CZy + Spy Eyy» 
1 
cwo = aa sez + gezz — €7 85, Cay — (97 — 9) Soy Czy] 
1 
cyo = Fag [29°cz + (g9? + 9” \ezz — cy — (9? — g”)cza] ; (A.3) 


where ôm can only be induced by custodial symmetry breaking effects and is set to 
zero in our framework. For the Yukawa couplings, we focus on the ones of t, c, b, T, y, 
parameterized as 
h 2 
Lnpp=—- J myp(l+ oys)frfr thc. (A.4) 
f=t,c,b,T, u 
The possible flavor violating Yukawa couplings from new physics are not considered. The 
anomalous triple gauge couplings (aTGCs) are parameterized as 
Lige = 1986 Al(W “Wy - W! *W iw) 
+ ig(1 + 597 )co, ZW" Wà, — WW) 
+ ig [(1+ 66z)CoyZ + (1+ dK) S64 AY] W; W3 


(A.5) 


v ph? 


1g v v — 
=F ma 0200 2 + Ay Sy A" \W, pwr 


where Vaw = ô Vy — V, for V = W*, Z, A. Gauge invariance further imposes the 
relations kz = ôgı,z — ty OK and Az = Ay. This leaves three independent aTGC 
parameters, which we choose to be 09),z, 0K and Az. Two of them, ôgı,z and ôk}, are 
related to the Higgs parameters and can be written as 


1 
091.2 = 2(g? — gi) -9(9° + g”)ez = gg + 9” )ezz + eg ox + g’ lg? — g’)cza| ; 


2 2 2 12 
g e g g 
OKy = — (e 2 g + er 24 git E ezz) a (A.6) 


To summarize, in our framework the contribution from D6 operators to the Lagrangian 
in Eq. (A.1) can be parametrized by the following 12 parameters: 


ÔCZ, CZZ, CZO, Cyy, CZy, Cgg, ÓYt, Ye, Yb, Ves. ÔYu, AZ. (A.T) 


Also following Ref. [3, 28, 24], we consider the EFT contribution to the hyy and hZy 
vertices at the tree level, in which case the only EFT parameter that contributes to the 
decay rate of h > yy (h > Zy) is cy, (czy). For the decay h —> gg, we include, in 
addition to Cyg, the contributions of dy, and dy, which enter the hgg vertex by modifying 
the Yukawa couplings in the fermion loops. It is also convenient to normalize Cy, Czy 
and Cgy with respect to the SM 1-loop contributions. We follow Ref. [3] and define the 
following parameters 


T as Vz = 
peat <1 — 26y, pat ~ l> 2¢zy, (A.8) 
yy Zy 
and r 
= ~ 14 200 ~ 14+ 2G +2.10 dy, — 0.10 dy, (A.9) 


SM 
Dga 
where Cy, Czy and Cg, are related to the original parameters by 


P m Cy Pp m o Zy P m _.. yg A 10 
Cr = g3 x 10-2” CZy = FO x 10-2” Ca = 83x 10-3" (a10) 


Furthermore, without measuring the tth process at high energies (/s = 500 GeV) or at the 
LHC, the parameters c,, and dy; can not be independently constrained. Since we focus on 
the 240-250 GeV run at lepton colliders, we replace cg, and dy, by Cc. in Eq. (A.9) which 
parametrize the total contribution to the hgg vertex. The number of parameters is thus 
reduced to 11, and the parameters are 


0cz, CZzZ, CZO, Cory 5 CZy, a. OYe , OY , OYr , OU Az, (A.11) 


which are used in our EFT global analysis in Section 3. 


B EFT expressions for the ete” — vvh cross sec- 
tions 


We obtain the cross section of ete” — vvh as a function of the EFT parameters by 
generating events using Madgraphd [14] with the BSMC package [29, 30]. The events are 
showered in Pythia [15] and passed to Delphes [16] with the ILD card for detector 
simulations, after which the selection cuts in Section 2 are applied. The interference 
between hZ and WW fusion are also included. The results for CEPC 240 GeV with 
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unpolarized beams are listed as follows. For the total rate, we have 
unpolarized 


ad = 1 +1.7őcz + 1.3czz + 2.9cz 


vvh |240 GeV 


+ 0.051 cz, + 0.146 


o 
+ 0.23 dew — 0.0026 cww — 0.065 cwo . (B.1) 


Here we do not impose the gauge invariance condition from Eq. (A.3) in order to show 
the different dependences on the Z and W parameters. For the binned differential distri- 
butions of Mrecoil (Msoi), the numerical coefficients in Eq. (B.1) are replaced by the ones 


in Table 4 (Table 5). We have also checked that the statistical uncertainties from sim- 


CEPC 240 GeV (with unpolarized beams) Mrecoil 
bin index [GeV] 

75 80 85 90 95 100 105 110 115 130 

osm [fb] 0.15 0.18 0.38 0.78 1.2 1.3 1.1 0.74 0.47 0.34 

ôcz 0.97 1.4 1.6 1.7 1.8 1.9 1.9 1.9 1.8 1.9 

CZZ 0.50 0.95 1.1 1.3 1.3 1.4 1.4 1.4 1.4 1.4 

CZ 1.5 2.2 2.6 2.9 3.0 3.0 3.2 3.1 3.1 3.0 
CZy 0.021 0.035 0.044 0.051 0.052 0.054 0.056 0.056 0.055 0.055 

Cy 0.075 0.11 0.13 0.14 0.15 0.15 0.15 0.15 0.15 0.15 
ôcw 0.93 0.62 0.37 0.22 0.15 0.13 0.14 0.17 0.20 0.29 
Cww -0.011 -0.0066 -0.0038 -0.0023 -0.0016 -0.0013 -0.0013 -0.0016 -0.0019 -0.0023 
cw -0.30 -0.18 -0.11 -0.060 -0.038 -0.032 -0.033 -0.036 -0.037 -0.049 


Table 4: The coefficients of the EFT parameters for the expression of o/agy of each bin 
of the Mrecoil distribution. The upper bound of each bin is listed in the first row. The 
first bin include all the events below 75 GeV. 


CEPC 240 GeV (with unpolarized beams) ar 
bin index [GeV] 

75 80 85 90 95 100 105 115 

osm [fb] 0.22 0.24 0.59 1.7 2.4 0.99 0.32 0.11 

ôcz 0.95 1.4 17 1.8 1.9 2.0 1.8 1.4 

CZZ 0.54 0.99 1.2 1.3 1.4 1.4 1.4 1.0 

CZ 1.6 2.3 2.6 3.0 3.2 3.3 3.1 2.0 
CZy 0.021 0.035 0.045 0.053 0.056 0.059 0.054 0.044 

Cy 0.075 0.11 0.13 0.15 0.16 0.16 0.15 0.11 

ôcw 0.92 0.61 0.33 0.14 0.075 0.11 0.33 0.99 
cww -0.0095 -0.0062 -0.0034 -0.0014 -0.00082 -0.0012 -0.0025 -0.0075 
cw -0.28 -0.17 -0.092 -0.037 -0.019 -0.025 -0.056 -0.13 


Table 5: Same as Table 4 but for the moi distribution. 


ulation are under control.” We then impose the gauge invariance condition in Eq. (A.3) 
and construct \? of the ete~ — výh measurement, assuming the events follow a poisson 
distribution. For the total rate, we have 
2 
N3, (1 2e) 
2 Oph 


aa Nsig + Nog í 


(B.2) 


TThere are nevertheless some small fluctuations in our result. For instance, the coefficients of ôcz and 
dcw should always add up to two. In most bins, the sum is controlled in the range 1.9-2.1. After imposing 
dcz = dcw we simply fix its coefficient to 2. We do not expect the fluctuations in other coefficients to 
significantly change our results. 
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where Ngig and Np, are the number of signal and backgrounds after cuts, normalized to 
5ab ' for CEPC. For the binned distributions, we use Eq. (B.2) to construct the x? of 
each bin where Ngig and Npg are the number of signal and backgrounds in the bin. We 
then sum over the y? of all the bins, assuming no correlation among them. The x? is 
then combined with the ones of other measurements for the global analysis in Section 3. 
We refer the readers to Ref. [3] for a complete set of expressions for the other relevant 
observables. 
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